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The use of the vacuum lepton pair production process (v ^ ve~ e'^), a, viable reaction for superlu- 
minal neutrinos, to put constraints on Lorentz violations requires a dynamical framework. Different 
choices of dynamical matrix elements and modified dispersion relations for neutrinos, leading to 
numerical factors differing by one order of magnitude in the results for the pair production decay 
width, are used to show the uncertainties on these constraints. 
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^ I. INTRODUCTION 

o 

Cn Lorentz invariance is one of the pillars of the current theories of fundamental physics. It implements the physical 

O , equivalence of inertial observers in absence of gravity, together with the presence of an invariant scale, the speed 

O of light in vacuum. In particular, Lorentz symmetry is an important ingredient of the Standard Model (SM) of 

C/j elementary particles. However, the logical consistency of the latter does not rely on the former, contrary to what 

QQ happens to unitarity or the absence of gauge anomalies which are necessary requirements for a quantum theory like 

^Nj the SM. This permits investigating departures from Lorentz invariance in the framework of a quantum field theory 

(QFT), and therefore, testing Lorentz invariance as a symmetry of the SM. Although much progress has been made in 

this direction, the non-compactness nature of the Lorentz group, the fact that the boost parameter can be arbitrarily 

O , large, implies that the unexplored region of this group will always be infinitely great, unlike the rotation group 

I which is compact and can be checked for every angle. Apart from the own interest of testing such a fundamental 

symmetry, there are motivations to study deviations from Lorentz invariance coming from both the theoretical and 

the phenomenological sides. On the one hand, some indications coming mostly from theoretical investigations of the 

quantum-gravity problem have suggested that Lorentz symmetry could be either violated or deformed at high enough 

energies (see [IHZj)- On the other hand, some observations might be interpreted to be incompatible with Lorentz 

\^ invariance, such as a possible violation of the GZK cut-off [8,, anomalies in the propagation of very high energy 

jj^ gamma rays in intergalactic space [TU] , or, more recently, the observation of apparent superluminal neutrinos by 

QQ the OPERA collaboration [TT following a proposal [l2] to look for signals of Lorentz invariance violation in neutrino 

ly-j propagation. 

^S| In this work, we present an analysis of the consequences of such departures from exact Lorentz invariance in 

-^ the decay of superluminal neutrinos. In particular, we shall focus on charged-lepton pair emission by superluminal 
f~^ neutrinos {vi — >■ Vi e" e^) which is a forbidden process in a Poincare invariant theory because of the charged 
^Nl lepton's masses. The idea that neutrinos can travel faster than light is not recent, and it actually goes back to the 

T-H eighties (see [13] and [HI). The difficulty in the detection of neutrinos (owing to the weakness of their interactions 

►^ with matter) entails some of their features remain still unknown, making neutrinos somehow mysterious particles. 
Since neutrino oscillations suggest that they are not massless particles, it is natural to expect the existence of right- 
handed neutrinos which are transparent to the electroweak and strong interactions of the SM (they are singlets 
under C/y(l) x SUl{2) x SUc{3) transformations). This special property has been used in the large-extra-dimensions 
scenario to speculate about the possibility that neutrinos could not be confined to our three-brane, but they are indeed 
exploring the extra dimensions [15^ . Because of this, if the departures from Lorentz invariance had their origin in 
these additional dimensions, it would be expected that neutrinos were the unique particles capable of directly feeling 
such deviations from exact Lorentz symmetry (perhaps together with an hypothetical graviton which would also be 
transparent to the gauge interactions of the SM) . The particular quantum numbers of neutrinos in the SM make also 
possible to write a gauge invariant but Lorentz violating Lagrangian which gives rise to effects of Lorentz invariance 
violation which are different for neutrinos than for the rest of the SM particles. Other works have considered neutrinos 
as the only candidates for Lorentz violating particles in an effective field theory (EFT) framework because in this 
context, and in presence of interactions, different limiting speeds of different particles at high energies are driven by 
the renormalization group flow towards a universal speed at low energies (although strong assumptions have to be 
made in order that the flow be sufficiently fast), the vacuum speed of light (see [H]). The weaker the interaction, 
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the weaker the running and therefore the difference of hmiting speeds is the greatest when the interactions are the 
weakest. This suggests that the hmiting speed of neutrinos at low energies may be different from the speed of the 
rest of particles when Lorentz invariance has been broken at high energies. 

In this paper, we shall parameterize the deviations from exact Lorentz invariance by a modified dispersion relation 
for free neutrinos. The consequence of such modification is twofold: on the one hand, there is a new dependence of 
the velocity on the energy which makes the time of fiight different from the one of a special-relativistic theory. On the 
other hand, the kinematics of particle processes involving neutrinos changes. In fact the dynamics of these processes 
is also changed because of the different choices of a dynamical matrix element compatible with the modified dispersion 
relation one can make. Previous works have focused on getting qualitative or quantitative results for the decay width 
Vi — )■ Vi e~ e"*" or Vi — )> Ui Vj i7j |17H21| . without discussing the different alternatives one has for the matrix element, or 
even without specifying the matrix element they are using in their calculation (however, see |22j). The main purpose 
of the present work is to show how this can affect the results for this process. One could think that a criterion to 
select a matrix element would be that it could be deduced from a local EFT. However, there exist limitations to 
this framework. First, one should restrict dispersion relations to analytical expressions so that a momentum power 
expansion be possible in the EFT. Second, the most likely origin of such deviations from Poincare symmetry arises in 
attempts to reconcile general relativity and quantum theory, through residual effects (that is, effects that are present 
when the classical curvature of space-time can be neglected) in the structure of space-time or/and momentum space 
which modify its classical, Minkowskian nature [231 [21]. If this were the case, it would not be surprising that this 
led to non-local effects at low energies since the semi-classical physics of black holes tells us that the fundamental 
degrees of freedom of gravity cannot be described by a local QFT, regardless of whether the latter has a non-trivial 
fixed point or not, because the densities of states of both theories do not match [321 HHj- Furthermore, doubts on the 
validity of an EFT description of the low energy limit of quantum gravity come from the difficulties to incorporate 
the necessary cancellations of contributions to the vacuum energy (cosmological constant problem) . 

Throughout this work, we shall assume: i) Rotational invariance is preserved, ii) Energy and momentum are 
conserved in the conventional, additive way. iii) The relevant propagation speed of the neutrinos is the group velocity 
of their wave packets, iv) Charged leptons are (or can be approximated by) special-relativistic particles, v) Indirect 
bounds on neutrino masses based on neutrino oscillations and cosmological observations are still valid, in such a way 
that neutrinos masses can be neglected in the process. 

This work has the following structure. In the next section (In]) we shall show the general procedure to obtain the 
decay widths, and we will compute them for four different choices of the matrix element and a general dispersion 
relation. Section |III| is devoted to discuss the properties of the previous matrix elements and their physical inter- 



pretation. In Section IV we shall specify dispersion relations to give definite analytical results for the decay widths 
and the rates of energy loss due to the pair emission. We shall also provide numerical results for these expressions. 
Section W\ is dedicated to the study of the time of flight of neutrinos in the presence of the aforementioned loss of 



energy. The following section (VI) is devoted to analyze the consistency of the super luminal speed interpretation of 



the observations reported by the OPERA collaboration. Our conclusions and remarks will close this work in Sec. VII 



II. PAIR PRODUCTION DECAY WIDTH 

Let us consider the process v(j)) — >■ v{pi) e~{p2) e^{p7i) induced by the production of a virtual Z^ and subsequent 
decay into an e~e+ pair. This is a kinematically forbidden process in special relativity (SR) which becomes allowed 
for superluminal neutrinos above a certain energy threshold. We are going to consider a very high energy neutrino so 
that all three particles in the final state are contained within a small cone around the direction of the momentum of 
the neutrino in the initial state. The decay width of the process is given by 
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where (|A^P) is the squared ampHtude averaged over initial spin states and summed over final spin states, 



{\M\')^A^^P,P,)P2^P,,, (2) 

and we are using the factorization of the matrix element into a factor depending on the neutrino variables and a 
factor depending on the electron-positron momenta. The second factor is the standard SR Dirac trace which gives the 
dependence on the four-momenta P2 and p^. The coefficients A^^'^ depending on the neutrino momenta p, pi contain 
all the corrections to SR in the dynamical matrix element. 



The decay width can be written as 
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r= ^ / 7,l^rj:t^^A^''{P^P^)B,AP~Pi) (3) 



where 



In the approximation where one neglects any modification to SR kinematics for the electron and positron, 

S^.(fc) - Bi{k^)7^^,e + B2{e)k^K , (5) 



where k^ — k^ — k'^, and 



1 / 4m 



2\3/2 



i..(fc^)^^(2 + ^)(l-5^)^^(.^-4™^).(M, (6) 

where 9 is the Heaviside step function. On the other hand, one has 

k' - ip-Pi)' = {E~E,r - (\p\ - |pi|)2 -2|pl|pi|(l -cos^i) (7) 

where ^i is the angle between the neutrino momenta p, pi . The minimum value of k^ {Am^ ) corresponds to a maximum 
value of (1 — cos 01 ) 

n cn.O) _ iE-E,r-{\p\-\p,\r-4ml 

(l-cos0i)+- ^^^^^ . (8) 

If we parameterize the energy-momentum relation for neutrinos aaj 

E=\p\il + e{\p\)) (9) 

then we have 

{E~E,f - i\p\ - \p,\)' = [iE+\p\) {E, + Ipil)] [\p\e{\p\) - \pM\pi\)] (10) 

and one has a negative result for the maximum value of (1 — cos^i), i.e., pair production is not allowed, in the case 
of SR kinematics (e(|p|) = 0). The threshold momentum \pth\ for pair production verifies: 

\pth\'e{\pth\)^2ml (11) 

corresponding to the configuration in which the outgoing neutrino is at rest, and the electron and positron have 
the same momentum. We shall consider incoming neutrino momenta much greater than this threshold momentum, 
\p\ S> Ipthl, in such a way that one can safely neglect electron and positron masses. In addition, we shall assume that 
the deviation from SR kinematics is sufficiently small (|e(|p|)| ^ 1) so that |pi|(l — cos0i)+ ^ |p|. In this case, the 
change in the direction of the neutrino energy flux due to the emission of the e^e"*" pair is very small ("coUinear" 
approximation) . 

We shall also restrict all the discussion to neutrino momenta such that 

[{E + \p\) (E, + \p,\)] [\p\ei\p\) \pM\p^\)] « M| . (12) 

In this case we can use the point like approximation for the fermion interaction in the pair production process and 
then the momentum dependence in the coefficients A'^'^{p,pi), which parametrize all our ignorance on the dynamics. 



^ We are neglecting neutrino masses in all the discussion. 



comes from the modified neutrino Dirac trace. Assuming a linear dependence on the components of each momentum 
in the amplitude then one has a general form for the angular dependence 

A'"'(p,Pi)B^,(p - pi) = ^^j^ [(1 - 2s^)2 + {2slvf] [^o(xi) + i^i(xi)(l - cos^i) + F2{x,){l ~ cos^i)^] ^^^^ 

X ((1 - cos ei)+ - (1 - cos^i)) e{l-xi+ e(|pl) - xie{xM)) , 

when one uses the momentum fraction xi — \pi\/\'p\ for the neutrino after pair production, and where Gp is the Fermi 
constant and sw is the sine of the Weinberg angle. The determination of the coefficients F„(a;i) of the three angular 
terms requires a definite dynamical framework. The angular integral on the neutrino momenta can be madcPl and the 
leading contribution for the decay width is 

-•2 1^5 /-If n _ r-^ofl. ^2 n _ ^r^clfl. ^3 1 

(14) 
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, ., ^ . ,,l-cos6'i)i , ,(l-cos6ii)3 
Fo(xi){l - cos0i)+ + FiixiY- ^-^^ + F2{xiy- ^— ^ 



X ^ (e(|pl) - xi€{xx\p\)) Xi dxi 
when one consistently retains only the dominant contribution in an expansion in powers of e(|p|). 

A. Dynamical matrix element 

Lacking a well defined dynamical framework incorporating superluminal neutrinos, we are going to consider four 
different simple choices for the dynamical matrix clement as a way to illustrate the uncertainties in the evaluation of 
the pair production process. Naively, one would take the SR matrix element as a first approximation and one would 
expect that any other choice approaching the SR limit would lead to the same leading order result for the decay width 
but this is not the case. The reason is that in the SR limit and neglecting masses all four momenta p, pi, p-, p+ 
are light-like and, due to energy-momentum conservation, proportional to each other so that any scalar product (and 
then the dynamical matrix element) vanishes. Then the leading contribution to the dynamical matrix element comes 
from the first non- vanishing correction to the SR limit. 

1. First example 

An apparently natural choice for the dynamical matrix element corresponds tcrl 

A^-^(p,p,) = 16G|. [(1 - 2s2^)Xp'^ -I- (2s2^)VK] (15) 

which is the expression of the SR dynamical matrix element but with p° = |p|(l -I- e(|p|)) and pi — |pi|(l -t- e(|pi|)). 
In this case one has 

A^-^{p,p,)B^,{k)^^ [{l-2sl,f + {2sl.r] [2{p,-p)e +A{p,-k){p-k)]e{k')0{ko). (16) 

Using the products 

2(pi ■p)^2 [\p\\p,\ {e{\p\) + e(|pi|)) + \p\\p,\{l - cos^i)] 

e = 2 m - bll) i\p\em - \pM\pi\)) - \P\\P1\{1 - cos^i)] 

2(pi . fc) = 2 m - \pi\)\pM\pi\ + IPil i\M\p\) - WMlPiD) + IpIIpiKI - cos^i)] 

2ip . fc) = 2 m - \p,mei\p\ + \P\ {\M\P\) - \pM\Pi\)) - IpIIpiKI - cos^i)] (17) 

we get 

Fo{x,) = 4 [x^{l - x^fei\p\) 6(xi|p1) + 2:ei(1 - x^) (e(|pl) + e{x^\p\)) (e(|p1) - Xie{x,\p\)) + ^i {<\p\) " ^i<xi\p\)) 

F^{x,) = 4 [-xj {e{\p\) + e{x,\p\)) + 3xi(l - x,) {e{\p\) - x,e{x,\p\))] , 

F2{xi)^-8xl. (18) 



^ The angular integral in (1 — cos^i) goes from to (1 — cos9i)-f, , except for those outgoing (very close to zero) momenta for which 

2 < (1 — cos6i)+, in which case the angular integral goes from to 2 . 
^ We shall neglect masses in all the matrix elements. 



When the angular integral is done one finds 
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[{I ~ 2slry + i2sl^f] / dx,eiem-^M^i\p\)) 



4x1(1 - x,)hmH^i\p\) (e(|pl) - ^ie(^ibl)) 
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6:ri(l - x,f {e{\p\) + e{x,\p\)) {<\p\) - Xie{xi\p\)f + (4xi(l - x,) + -(1 - x^f) {e{\p\) ~ x,e{xi\p\)) 



(19) 



If one is not interested in the angular dependence of the differential decay width, there is a simpler way to get the 
result in (19). Going back to (16), one can use the variable fc^ instead of the angle 6i to express the products 



and 



2(P1 ■p)^p''+pl-k^ 2(pi . fc) = p2 _p2 _ ^2 



G 



2{p-k)^p'^ -pl + k^ 



A^'^ip,p,)B^,ik) = ^ [{1^2slf + {2sl^f] [ip^-plf + ip'+pl)k'~2{kY] e{k')eiko). 



(20) 



(21) 



The relation (l7| allows replacing the angular integral by an integration over k^, where the upper limit on fc^, corre- 
sponding to 01 = 0, is given by 



kl = 2\pWl-Xi)[e{\p\)^Xie{x, 



This leads to: 



Using 



^-i^t(i-2^-)^-^(^^-)^]2^r^-^(^^) 



p' = \p\'{l + 2e{\p\)) 



ip'~pifki + l{p'+pi){kir-l{kir 



p? = 1^112(1 + 26(1^11)) 



(22) 



(23) 



(24) 



one recovers the result ( 19 ) for the decay width 



The choice (15) for the neutrino dependent factor in the dynamical matrix element cannot always be derived from 



a quantum field theoretical calculation. In fact, by considering a generic dispersion relation (arbitrary choice of 
e(|p|)) one can have cases where Eq. (16) and the expression of the decay width take negative values indicating an 
inconsistency of the ansatz for the dynamical matrix element. 



2. Second example 



One could consider other alternatives to (15) for the matrix element. A very simple choice corresponds to consider 



a modified neutrino spinor satisfying a modified Dirac equation 

[7°i?(b1)-7-p(l + e(b1))]«(p)=0. 



(25) 



This modified Dirac equation implies a modified dispersion relation i5(|p|) = |p|[l + e(|pl)]. With this modified 
Dirac neutrino spinors the matrix element can be calculated as in SR and the result for A^^"^ is now 



i''''(p,Pi) = 16 G| [(1 - 2sl.fp^,p-' + (2s2^)Vp5'] 



(26) 



with 



/ = b1[l + e(b1)] P^p[l + e{\p\)]. (27) 

This neutrino factor can be derived from a perturbative field theory calculation by considering the SR vertex for the 



interaction and a modified free fermion action leading to the simple modification of the Dirac equation ( 25 ) 
With the choice ( 26 ) for the neutrino factor in the dynamical matrix element one has 



i^'^(p,Pi)S^,(fc) = 16G|,[(l-2s^)2 + (2s2^,)2] [{p,.p)k^B,ik^) + (pi.fc)(p.fc)B2(fc2)] 

Q2 (28) 
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[(1 - 2sl^f + (2s2^)2] [2(pi • p)fc2 + 4(pi ■k){p- k)] 9{k^) 9iko) . 



In this case we have the products 

2pi-p = 2|p1|pi|(l-cos0i) 

e = 2 [{\p\ - \p,\) {wm - \pMm)) - biipii(i - cos^i)] 

2(^1 ■ k) = 2|pi| iWrn) - \pM\pi\)) + 2|p1|pi|(1 - cos^i) 
2(p. fc) = 2|p1 (\M\p\) - \pM\pi\)) ~ 2|p1Ipi|(1 - cos^i) (29) 

and 

Fo(a^i)=4a;i(e(|pl)-xie(a;i|pl)f , 

Fi(xi) = 8x1(1 - xi) (e(|pl) - a;ie(a;i|:?^)) , 

F2(2;i) = -8x?. (30) 

The decay width after the angular integration takes a different and simpler form than the decay width result of the 



first example ( 19 ) 



f = ^1^ [(1 - 2s^wf + i'^^wf] \ I' dx, e (6(b1) - x,e{xM)) (1 " ^?) {<\P\) ~ x,e{xM)f ■ (31) 

Note that the result for the dynamical matrix element and then for the decay width is positive definite independently 
of the choice of the modified dispersion relation, so that the potential inconsistencies of the use of ( 15 1 for the neutrino 
factor are not present in this second example. This is a direct consequence of a field theoretical derivation of the 
dynamical matrix element. 

3. Third example 

One has a third simple choice for the dynamical matrix element with 

A^'^ip.p^) - 16G| [(1 - 2sl,fp'tp'^ + (2s2^)Vp5'] (32) 

with jP — \p\ and p — p which is just the SR dynamical matrix elementrj All one has to do is to replace in the first 
calculation everywhere p by p. Then one has 

A^^%p,Pi)B^,{k)^l&Gl[{l-2s''^f + {2slrf] [{pi-p)k'B^{e) + {p^ ■ k) {p ■ k)B2{k')\ 

n2 (33) 

= Y^ [(1 - 2sl,f + {2sl,f] [2(p-i . p)e + 4(pi ■k){p- k)] e{k') e{ko) . 

Due to the proportionality p^ = P^(l + e(|p|)) it is obvious that the decay width in this third example coincides with 
the decay width in the second example up to corrections of order e"*. 

4- Fourth example 

In the particular case of a modified dispersion relation for the neutrino with a momentum independent choice for 
e, which corresponds to a momentum independent speed, it is possible to consider a modification of the interaction 
vertex fixed by gauge invariance from the modified free fermion action. If we consider a Lagrangian 

C = DLij"Doi^L-i{l + Y)i^Ll-DiyL (34) 

where Dq, D are covariant derivatives, then one has the simplest generalization of the relativistic Lagrangian for 
gauge interactions of massless fermions with a constant velocity v = 1 + ?7o/2. In order to derive the dynamical matrix 
element all one has to do is to replace the Dirac 7 matrices j^ in the neutrino tensor by the modified 7 matrices 7'' 



Note that in the first choice of the dynamical matrix element one takes into account the modified expression for the energy at the level 
of the matrix element. 



where 7° = 7° and 7 = (1 + 770/2)7. This replacement incorporates the modification in the Dirac equation for the 
modified Dirac spinors and the modification in the gauge interaction. Using 

{r.ry^^r" (35) 

where t)^" is the modified Minkowski metric fy™ = 77'^'' = 1, 77" = (1 + 770)77" = — (1 + 770) one can show that 
The coefficient of Bfj_i,{k) in the dynamical matrix element can be read in this case from 

PipP. [v'^r^ - f]'"ff^ + ^"^rl [k^i - vcpv'-" + s:ls'^) (37) 

which is the symmetrized coefficient of P2f_iP3i> in the product of the two traces. By contracting indices one finds 



275^755' + 27)5'75" + (4 + 3770) (75i • p)??^" - 2(75i • 75)^7^" - 2(7Ji • p)??^" , 
to be compared with the SR coefficient 

PipP. W"V"" - V'^r^"^ + ^"^'^1 {K^^'p - Vcfsv'^" + S:;,sfj = 2p>r + 2p^,p'' . 
Then the dynamical matrix element will be 

A'^^{p,p,)B^,{k) = ^ [(1 - 2^2^)2 + {2sl.r] [-4{p, ■ p)e + 4(7)1 • k){p ■ k) 

+ 2 (4 + 3770) (731 • p)fc' - 2(7)1 • p){k ■ fc)] e{e) e{ko) , 

with 

p = {\m + ^o/2),p*(i + %)) P = {\m + ^0/2), 7)) 

7)1 - (|7?i|(l + ??o/2),7?i(l + r/o)) Px - (|pi|(l + %/2),7?i) 

fc = ((N-|Pi|)(l + '?o/2),(7?-pi)(l + r?o)) fc-((|pl-|p'i|)(l + r7o/2),(7?-7?i)). 

At first order in an expansion in powers of the corrections to SR one has 

2(7)1 -75) = 2|p1|7?i|(1 - cos^i) - 2|pl|7/i|7;o 

e - -2|7Jl|7?i|(l - cos^i) + (IpI - |7?i|)'r;o 
2(7)1 • fc)-2|73l|7?i 1(1- cos 0i) 
2(7).fc)--2|77l|77i|(l-cos0i) 
2(7)1 • p) = 2|p1|7?i 1(1- cos 0i) 
(fc-fc) = -2|7jl|77i|(l-cos^i) 



and then 



h{x\) =4xi(l-xi)2772, 
A(a:i)= [4xi(l-xi)2-8x2]77o, 



After integration on the angle Q\ , we get 
f 



G|bT 



[(1 - 2sl.f + {2sl.f] 770^ / dxi 

Jo 



x,{l - x,)^ + -{1 - x,)^ 
6 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



which is a third candidate for the pair production decay width, this one limited to the case of a superluminal neutrino 
with a momentum independent speed v = 1 + 770/2. 



III. CHOICES OF DYNAMICAL MATRIX ELEMENT AND MODIFIED DISPERSION RELATION 

FROM A FIELD THEORETICAL PERSPECTIVE 

From a field theoretical perspective one should start by considering a generalization of the Lagrangian of the SM 
containing the neutrino field. From the free part (quadratic in the neutrino field) one could read the generalization of 
the Dirac equation and the corresponding modified dispersion relation (i.e. the function e(|p|)). By also considering 
the interaction term (product of two neutrino fields and the Z-boson field) one could derive the expression for the 
dynamical matrix element ((|A^p)). 

Since one is considering a generalization of a relativistic gauge theory the most natural way to implement a mod- 
ified dispersion relation in an extended Lagrangian is to add a new term in the free fermion action that fixes the 
modification in the dispersion relation and to consider the same dynamical gauge principle (replacing partial by co- 
variant derivatives) that fixes the gauge interaction in the relativistic limit. The fourth example for the dynamical 
matrix element considered in the previous section is just the simplest choice along these lines with just one spatial 
derivative in the free part. One could consider generalizations with higher spatial derivative terms that would lead 
to new interaction terms and new contributions to the dynamical matrix element. A generalization along these lines 
is restricted to the effective field theory framework (derivative expansion) so that the choice of a modified dispersion 
relation is very limited and the study of implications of (or constraints on) such a generalization is easier than in 
other cases. To be precise, the gauge invariance that has been implemented in the Lorentz violating generalization of 



the SM Lagrangian in Eq. (34) is a C/(l) symmetry instead of the complete SUl{2) x C/y(l) gauge invariance of the 
Lorentz invariant contribution. Should we had considered an SUl{'2) x tV(l) gauge invariant extended Lagrangian, 
Lorentz violations on different particles would be strongly restricted. In particular since the neutrino field and the 
left-handed charged lepton field are in a gauge doublet then Lorentz violations in the neutrino sector automatically 
would have an analogue in the charged lepton sector. The absence of observations of such Lorentz violating effects 
for electrons implies strong restrictions on posible effects due to Lorentz violations in neutrino physics including the 
possible energy loss of superluminal neutrinos that we are discussing in detail in this work. 

An alternative way to implement a modified dispersion relation is to consider an extension of the relativistic La- 
grangian independent of the gauge fields constructed from gauge invariant products of matter fieldsj^ This possibility 
is restricted by the matter field content of the theory. In the case of the SM Lagrangian the fact that the (conjugate 
of the) doublet scalar field $ and the left-handed doublet lepton field L 

have the same gauge transformations allow to consider an extended Lagrangian quadratic in the gauge invariant 
product of these two field doublets 

CLivi^^L) = -^ (Zl.) i7 • V e(|zV|) (l>ti) . (46) 

In the approximation where one neglects the fluctuations in the scalar field this extended Lagrangian reduces to a 
quadratic Lagrangian in the neutrino left-handed field. This opens the possibility to consider Lorentz violating effects 
in the neutrino sector with no analogue for other particles so that restrictions from the absence of observations of 
such effects in other systems do not translate directly into restrictions on possible effects of Lorentz violations in the 
neutrino sector. This provides us with an example where neutrino physics is a special window to explore departures 
from SR. As far as one does not require a perturbative treatment of the Lorentz violating interactions in the lepton- 
scalar sector there is not any restriction on the quadratic extended Lagrangian so that one can consider arbitrary 
choices for the modified neutrino dispersion relation (arbitrary choice for e(|p|)) going beyond a derivative expansion 



(momentum power expansion) . The second example ( 26 1 for the matrix element of pair production by superluminal 
neutrinos can be seen as a result within this framework. 

A third alternative to the generalization of a relativistic gauge theory would correspond to assume that the local 
gauge symmetry is a property of the relativistic limit. In this case one could consider an extended Lagrangian with 



no restrictions from local gauge invariance. The second example (26) could also be seen as a tree level approximation 
within this framework which can not be distinguished from the second alternative unless one goes beyond this ap- 
proximation. An argument in favor of the realization of the second example for the dynamical matrix element within 
the previous second alternative for a generalized relativistic gauge theory is that it is not clear how the introduction 



^ The consistency of this field theoretical framework in the presence of radiative corrections can be questioned [27| and deserves further 
study. 
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of Lorentz violating terms in a Lagrangian provides a way to escape to the inconsistencies of a gauge non-invariant 
relativistic field theory with vector fields. 



The first example (151 for the dynamical matrix element can not be derived from a field theory perturbative 
calculation. It therefore requires to consider a generalization of the relativistic gauge theory that goes beyond the 
field theory framework and then there is no reason to consider restrictions on the choice of modified dispersion relation. 



For the second example ( 26 1 one can consider a momentum power expansion for the modified dispersion relation 
(e(|p|) if one assumes the validity of the effective field theory framework for the study of Lorentz violating effects or 
a more general momentum dependence if one assumes that one has to go beyond the effective field theory expansion 
when one goes beyond the special relativistic limit. 

The results for the third example, although can not be derived from a field theory calculation, are equivalent to 
those of the second example at leading order in the deviations from SR and then do not require any independent 
discussion until one goes beyond the leading order effects. 

Finally in the fourth example for the dynamical matrix element there is a definite modified dispersion relation 
corresponding to a momentum independent speed (in the massless limit). We could go beyond this case by including 



gauge invariant higher derivative terms in the Lagrangian ( 34 ) 



IV. ENERGY LOSS OF SUPERLUMINAL NEUTRINOS 

In this section we evaluate the width for pair production, the rate of energy loss, and the energy of a superluminal 
neutrino after propagation over a given distance for different choices for the dynamical matrix element and the 
modified dispersion relation. These results are the starting point of an analysis of the observable consequences of 
having superluminal neutrinos with the uncertainties due to the lack of knowledge of the details of the possible origin 
of the Lorentz invariance violation in neutrino physics. 

A. Decay width 



We have three candidates (191, (31 1, (44 1 for the decay width as a functional of the modified dispersion relation 
corresponding to different choices for the dynamical matrix element of the pair production process. 

If we consider the simplest choice for the modified dispersion relation e{\p\) = ?yo/2 (momentum independent speed) 
with 770 > (for negative values of 770 the decay width is zero) then the threshold is \pth \ = \/^n\\jr\Q and the decay 
widths are: 



ro = ^r^ [(1 - ^^'wf + i^-sl^r] eo ^ (48) 



To 
with 



192 7r3 



[il^2s'^r + i2sl^r]io{^y (49) 



.8 e _ 34 P _ 16 ,^„x 

^"~7 ^°-105 ^"-35- ^^^> 



The result for the decay width (47) corresponding to the first choice for the dynamical matrix element (SR matrix 
element with the replacement of the SR energy by their modified expression in terms of the momentum) reproduces 
the result used in Ref. |18j to argue against the consistency of the recent result of OPERA [TT] for a superluminal 



neutrino velocity. The other two results in Eq. ( 50 ) also reproduce the two cases studied in Ref. [22] , which in fact can 
be seen to correspond to the second and fourth matrix elements of the previous section for a momentum independent 
velocity. Our more direct computation will however allow us to consider other dispersion relations beyond the constant 
speed case. 



We can see from ( 50 1 that a change in the choice of the dynamical matrix element produces an additional overall 
factor of (17/60) in the decay width if one takes the dynamical matrix element corresponding to a perturbative 
field theory calculation with the simplest free action implementing a modified Dirac equation and the SR interaction 
(second example of the previous section) or an overall factor of (2/5) when one uses an interaction fixed by the 
dynamical gauge invariance principle (fourth example). 
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The next choice we consider for the modification in the dispersion relation is e(|p|) — |p|"/A" (with n and A positive 
numbers) which corresponds to a Lorentz violating free term in the Lagrangian with n spatial derivatives when n is 
a natural number. The energy scale A is the UV scale that fixes the domain of validity of the effective field theory 



energy power expansion. In this case, the threshold is given by \p-. 



'th\ 



|n+2 _ 



= 2 m^ A" and the widths are: 



G|b1 



192 



with 



Cn=2- 



"" 192^3 



12(n + 6) 



^ [(1 - 2s%,f + (2.^^)^] (M^ ^„ 
[(l-2.^f + (2.2^f 



Cn 



(51) 
(52) 



12(?i + l) 



4(3n + 2) 



(n + 3)(n + 4)(n + 5) (n + 2)(2n + 3)(2n + 5) (3n + 4)(3n + 5)(3n + 7) 

12 6 4 

(n + 2)(n + 5) ^ (n + 3)(2n + 3) " (3n + 4)(3n + 7) ' 



(53) 



The main difference between (511-^52^ and the results (47)-(48) for the simplest choice for the dispersion relation is 



the power exponent in the momentum dependence of the decay width. This makes the effect of the production of 
pairs in the propagation of neutrinos to increase much faster when the energy increases. 

As a third choice for the modification of the dispersion relation at high energies (|p| > po)j we use e(|p|) = A"/|p|" 
witq^O < a < 1 and A an additional energy scale required by dimensional arguments. This corresponds to a nonlocal 
free field theory action trying to illustrate a Lorentz violation in the neutrino sector that goes beyond the effective 



field theory framework. In this case, the threshold is given by |p^ 



'th\ 



2rril/X°' and the widths (for \p\ ^ |po|) are: 



192 7r3 

Gl\pf 
192 7r3 



[(1 
[(1 



with 



La 



12(6 -a) 



(3-a)(4-a)(5-a) 
12 



2sl^f + i2sl.r] 
2sl.r + (2.^)^] 

12(1 -a) 
(2-a)(3-2a)(5-2a) 
6 4 



3a 



^-c 



3q 



e-c 



(54) 
(55) 



4(2 - 3a) 



(4-3a)(5-3a)(7-3a) 



(2-a)(5-a) (3-a)(3-2a) (4-3a)(7-3a) 



(56) 



In this case the increase of the decay width of pair production with the energy is slower than in the case of a 
momentum independent velocity of propagation. In addition, note that in the expression of the decay rate there are 
energy- independent, dimensionless factors (i.e., pure numbers) that depend on the particular form of the dispersion 
relation and which can differ by one order of magnitude. For example. 



6 



^-3/4 



18.8. 



(57) 



As a last illustrative example, we shall consider a case in which the dispersion relation is non-analytic and it is 



jiven by e(|p|) 



-A/\p\ 



In this case, that also goes beyond EFT, the decay width is extremely sensitive to changes 



in the momentum. When \p\ <^ A, the dominant contribution to the decay width for the second matrix element is: 






(58) 



While this expression is small in the region in which is valid, it illustrates the possibility of having a remarkably strong 
dependence on momentum. 



If o > 1 then one has subluminal neutrinos and no pair production. 
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B. Rate of energy loss 

An approximation to the effect of the production of e+e^ pairs on the propagation of neutrinos can be obtained 
from the rate of energy loss 



m 

dt 



-\p\J^ d.i(l-.i)^ 



In the case of e(|p|) — 770/2 we have 

d\p\ 



2 1^6 



GI\P\ 



dt 192 7r3 

dt 

d\p\ ^ G|bt 

dt 192 7r3 



[(1-2. 



,2 \2 
'W) 



{^■^'wf] eo (I)' 



192.3 [il'2sl.r + i2sl.r]i'o''^" 



[(1 - 2.^^)2 + (2.?,)2]e'o (I)' 



with 



^0 



25 

28 



fo 



11 

42 



for the decay width results in (19), (31), and (44), respectively. 



If one considers e(|p|) — \p\"J7i^ then one has 

d\p\ Gl\p\^ 



dt 

m 

dt 



192 7r3 

G%\P\' 
192 7r3 



[(1-2, 



,2 \2 



2 \21 



i^siy) 



fo 



14 



3n 



c; 



^'^"^^'[(1-2.2^)2 + (2.2^)^ 



3n 



with 



c = 



12(2n + 15) 



12(4n + 3) 



e. 



13 

10 " (^3)(n + 4)(n + 5)(n + 6) ' (2n + 3)(2n + 4)(2n + 5)(2n + 6) 

2(5n + 4) 10 

~ (3n + 4)(3n + 5)(3n + 6) ^ 3(3n + 7)(3n + 8) 

3 24(n + 4) 24(2n + 5) 



(59) 

(60) 
(61) 
(62) 

(63) 

(64) 
(65) 

(66) 



5 (n + 2)(n + 3)(ri + 5)(n + 6) (2n + 3)(2n + 4)(2n + 6)(2n + 7) 
8(3n + 6) 



(67) 



(3n + 4)(3n + 5)(3n + 7)(3n + 8) ' 
With the third choice of modified dispersion relation, e(|p|) = A"/|p|", the rate of energy loss is given by 



m 

dt 

m 

dt 



Gi\p\ 
192 



^ [(1 - 2sl,r + i2sl.f] (^ 



2 1,^6 



GI\P\ 
192 7r3 



[(1 - 2sl.r + i2sl.)'] 



> — a 



f- 



with 



e-c 



12(15 -2a) 



12(3 -4a) 



13 

10 " (3^a)(4-a)(5-a)(6-a) ' (3 - 2a)(4 - 2a)(5 - 2a)(6 - 2a) 

2(4 -5a) 10 

~ (4 - 3a)(5 - 3a)(6 - 3a) ^ 3(7 - 3a)(8 - 3a) 

~ _ 3 24(4 - a) 24(5 - 2a) 

^'-" ~ 5 ~ (2-a)(3-a)(5-a)(6-a) ^ (3 - 2a)(4 - 2a)(6 - 2a)(7 - 2a) 
8(6 - 3a) 



(68) 



(69) 



(70) 



(4 - 3a) (5 - 3a) (7 - 3a) (8 - 3a) 
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C. Final energy after propagation 



The previous results of the rate of energy loss can be used to get an estimate of the final energy Ef of a neutrino 
of energy Ei after propagating over a distance L. 
One has 



1 1 



1 



1 1 



1 



with 



Ej Ef El 



En 



Ej Ef El 



1 11 

^~^~Ei 



§^[ii-2sirH2s^^r],{^yir' 



(71) 



(72) 



in the case of a momentum independent velocity (e(|p|) = '7o/2)- The result for Eq (Eq) is obtained from (72 1 by the 
replacement of ^'g by ^g (^'o)- 

For a modification of the dispersion relation with e(|pj) — |p|"/A" one has 



1 



1 



p5+3n rnS+Sn 



E:^ 



5+3n 



pS+Sn p5+3n 



En 



■5+3n 



(73) 



with 



E„ 



^^^ [il ^ 2sl.r + i2si^f] (5 + 3n)-i^e; 



192 7r3 



A3" 



-l/(5+3n) 



(74) 



and a similar result for En replacing ^'„ by ^^. 
Finally, in the case e{\p\) — A"/|p|" one has 



1 



1 



7715— 3a p5 — 3a 



E 



5-3q 



IT'S — 3a 7^5— 3a 



E 



5-3a 



with 



(75) 



E-a = 



GlL 



-,-l/(5-3a) 



192 7r3 
and E_a with a factor S^'_^ instead of f'_ 



[{l-2slrf + {2sl,f] (5-3a)A3"Ca 



(76) 



D. Some numerical estimates 

We can take the inverse of the decay width of pair production as an estimate of the distance that neutrinos should 



propagate to have an appreciable loss of energy. From ( 48 1 we have 



~ 1 _ 4.3 X 10-'' /GeVV 



^ll 



km 



(77) 



and similar results for Fq ^ (f q ^) with an additional factor 17/60 (17/24) in the inverse of the decay width. In the 
case of a Lorentz violating free term with n spatial derivatives one has 



_ 1.7 X 10(5'^"-5) 

" - 1 



A 



GeV\ 



1019 GeVy \\p\ J 



( 



(5+3n) 



km 



(78) 



and in the case of a Lorentz violating correction which decreases at large momenta (55), 

1.7 X 10(27"-5) /eV^ '" 



TZ = 



^ 



A y [\p\ J 



km . 



(79) 



13 




hlEp 



E2/EeW 





















L^ 


1 pc (1" matrix element) 


0.6 


: 





L^ 
L^ 


I Mpc (1" matrix elemem) 
1 pc (2"*^ matrix elemem) 






0.5 


: 




L = 


1 Mpc (2"*^ matrix element) 




0.4 


: 








0.3 


: 






0.2 


: 




/ 


0.1 


- 






^^ 



0.001 



0.01 



0.1 



10 



A/Ep 



FIG. 1. Value of the terminal energy En (approximate final energy of a neutrino with initial energy higher than E„) for the cases 
n = 1 (left) and n = 2 (right) as a function of the high-energy scale A written in terms of the Planck energy Ep. The curves 
show the differences with respect to the matrix element used in the calculation and with different distances of propagation, 
typical of galactic or extragalactic neutrinos. 



These expressions give us an idea of the sensitivity to an energy-loss due to pair production in the propagation of 
neutrinos from observations of solar, reactor, accelerator or atmospheric neutrinos. 



The results for the energy scales (72), (74) and (761, 

1/5 



Eo 



1.6 X 10" 

3/5 



km 
17 



GeV 



En — 



E^r 



'1.7 X 10(57"-5)\ l/(5+3n) 

. (5 + 3n)4 J 

'l.7 X 10(27"-5)' 

(5 - ia)to. 



A 



1019 GcV 



3n/(5+3n) 



k^xl/(5+3n) 



l/(5-3a) 



eV\ 



3Q/(5-3a) 



km 



l/(5-3a) 



GeV 



(80) 
(81) 
(82) 



vifhich allow us to determine the final energy after a propagation over a distance L, give us also another estimate of 
the possible constraints that one can get on (or hints of) Lorentz violating corrections to neutrino physics from the 
observed high energy neutrino spectrum. 



examples of matrix elements one has 

1/5 



The results in (|80|)-(|82|) correspond to our second example (31) for the dynamical matrix element. For the other 

1/5 



Eq = Eq 






Eq = Eq 




En = E„ 




l/(5+3n) 



E. 



S-. 






l/(5-3a) 



For galactic and extragalactic neutrinos, it is convenient to change units and rewrite, e.g., Eq. (77) as 

5 



r, 







1.4 X 10^ 



-32 



/TeV\^ 



and 






1.4 X 10- 



'7o 



10^ GeV 



b1 ) 



Mpc. 



(83) 
(84) 

(85) 



As an illustrative example. Fig. Ill gives the value of En for the two most studied cases in the literature of quantum 
gravity phenomenology, n — \ and n = 2, in the case of galactic (propagation distance of the order of pc) and 
extragalactic (propagation distance of the order of Mpc) neutrinos. We also show in that figure the differences in the 
use of the first or second matrix elements in the calculation. 
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V. TIMES OF FLIGHT OF SUPERLUMINAL NEUTRINOS 



A modification of the neutrino dispersion relation produces an energy loss in the propagation of neutrinos related 
with the time of flight. Then one can try to look for neutrino observations where one can simultaneously determine 
the energy loss and time of flight in order to test the consistency of the observations with a given modification of 
the dispersion relation. The determination of the energy loss is limited by the uncertainties of our knowledge of the 
sources of very high energy neutrinos. As for the time of flight determination one has together with the uncertainties 
in the knowledge of the time of emission also the uncertainties due to the limited precision of time measurements. 

In the case of a momentum independent superluminal velocity it is straightforward to calculate the time of flight 
of a neutrino propagating over a distance L 



or 



Atn = 3.3 X lO'^ 



Atn = 1.3 X 10" 



(^)( 



1- — 



/ L 



VMpc 



)( 



1-^ 



(86) 



(87) 



A deviation, with respect to the SR expectation, in the time of flight of neutrinos of t he o rder of 10 ns in the 
propagation over a distance of the order of 10'^ km would correspond to 770 ~ 10~^. From (771 one can see that an 
inverse decay length of 10'^ km for such a value of 770 corresponds to a neutrino with \p\ ~ 40 GeV. This provides us 
with a quantitative comparison of the sensitivity of neutrino times of flight and spectrum observations to a departure 
from SR kinematics with a momentum independent speed for the neutrinos. Times of flight of extragalactic neutrinos 
are much more sensitive to deviations from SR with a deviation of the order of a few s for a neutrino propagating over 
a Mpc corresponding to much smaller values of 770 (% ^ 10~^^). An inverse decay width of a Mpc with this value of 
rjQ corresponds to a neutrino with \p\ '^ TeV. An extragalactic neutrino flux extending above these energies excludes 
then such a deviation from SR. 

In the case of more general modifications of the dispersion relation the time of fiight of neutrinos will be affected by 
the energy loss due to the production of e^e~ pairs. The momentum and then the velocity of the neutrino is changing 
in the propagation according to ( 65 ) or ( 68 ) . It is very easy to calculate in these two cases the neutrino time of flight 



over a distance L. In the case of e(lp|) = Ip|"/A" one has 



Ai„ = L- L{n+l) 



with 



(5 + 3n) 

(5 + 2n) 



m 



(5+3n) 



1- 



1- 



A 



El 

En 



(5+3n)l (5+2")/(5+3«)N 



In the case of e{\p\) = A"/|pj" one has 



At_„=i-i(l-a) ( -^j T_ 



with 



(5 -3a) fE^ 



(5 - 2a) V Ej 



(5-3q) 



3l 

E-a 



(5-3q)' 



(5-2a)/(5-3a)N 



(89) 



(90) 



(91) 



In the case oi E^r ■ L <$^ 192tt'^/GJp one has Ef ^ E^ {Ef <C £'-q) and r„ « 1 (t_q « 1), i.e., the time of flight 
corresponding to a uniform motion with a momentum dependent speed. 

VI. CONSISTENCY OF THE SUPERLUMINAL INTERPRETATION OF OPERA OBSERVATIONS 



Recently there has been a claim of an observation of neutrinos propagating with superluminal velocities in the 
CNGS beam from CERN to Gran Sasso ^1^ (OPERA experiment). It is then natural to try to accommodate these 
observations within the present discussion of the effect of Lorentz violations on the propagation of neutrinos. 



15 

Almost all the theoretical discussions related to OPERA consider a dispersion relation with e{\p\) = 770/2. This 
is due to the absence of a change in the measured times of flight in all the range of energies detected going from 
10 GeV up to 100 GeV. These results, if confirmed with more statistics, put strong constraints on any momentum 
dependence of the velocity of propagation of neutrinos at least in the range of energies covered by the OPERA time 
of flight measurements. But in fact the OPERA value for 770 = 4.7 x 10^^ [11. is in obvious confliclF] with time of 



flight limits from SN1987A [28H3T| as one can see from Eq. (87). This conflict requires to consider a deviation from 
a momentum independent choice for e at energies below those explored by OPERA including the range of energies 
of SN neutrinos (few MeV) . There is also a conflict of this large value of rjo with high energy neutrino observations 
owing to its implications on the propagation of high energy atmospheric neutrinos. This requires a strong suppression 
of the production of e+e" pairs at energies well above those explored by OPERA with an appropriate choice of the 
momentum dependence of e. This is also needed to avoid incompatibilities with neutrino production owing to the 
kinematic corrections in the pion decay induced by such large value of rjQ [32H34| . 

The simplest way to try to escape to these contradictions is to use a modified dispersion relation combining a 
function e{\p\) which increases when going from SN momenta to OPERA momenta, stays almost constant over the 
energy range of OPERA observations to be consistent with the very mild momentum dependence of the time of flight 
of neutrinos and then decreases if we go beyond OPERA momenta in order to escape to incompatibilities with the 
observed high energy neutrino spectra. However, there is still another conflict owing to the effect of the production 
of e~'"e~ pairs in the propagation of neutrinos from CERN to Gran Sasso. According to the previous section, one can 
see that, for the value of 770 required to reproduce OPERA results on time of flights, there is a very drastic correction 
on the neutrino energies in conflict with observations. 



For a momentum independent velocity of propagation in the OPERA energy range, taking Eq. (80) with 770 — 
4.7 X 10^^ and L — 731km as the distance traveled by neutrinos in the experiment, we get Eq — 17.0 GeV in the 
case of the second matrix element. In the case of the first matrix element, which is the one considered by Cohen and 
Glashow in Ref. [TS] , we obtain Eq = 13.3 GeV (in fact they get a slightly different result, Eq = 12.7 GeV, because they 
use 770 = 5.0 X 10~^, a value which was updated by the OPERA collaboration in November 2011 to 770 — 4.7 x 10"^, 
which is the one used in our calculation). For the third matrix element, Eq = 16.0 GeV. This means that together 
with the rather adhoc choice of the Lorentz violating corrections in the free Lagrangian which would produce the 
peculiar behavior of the e{\p\) function as described in the previous paragraph, one has to assume a suppression in 
the dynamical matrix element relevant for the calculation of the energy loss of neutrinos propagating from CERN to 



Gran Sasso. From Eq. ( 72 ) , we get 



^° " °-^^^ (imi) ''o"'^' ^""'^' ^""^ ■ ^^^^ 

Let us consider this relation for a generic matrix element. Since OPERA observes the arrival of neutrinos of energies 



higher than, let us take, Eq > 60 GeV, Eq. (92) gives an upper bound for ^'q-. 

Co < 4.8 X 10"^ (93) 

This means that in order to make the OPERA time-of-flight measurement compatible with the non-observation of pair 
production, one needs a suppression in the dynamical matrix element so that £,' ^ 0(10^^) instead of the 0(10^^) 



values of Eq. ( 63 ) obtained in the three simple examples considered in this work. 

One could wonder whether relaxing the OPERA observation of a constant velocity in the 10-100 GeV energy range 
could make the strong constraints from neutrino decay compatible with the value 770 = 4.7 x 10^^ measured at an 
energy oi E = 17 GeV (mean energy in the OPERA experiment) Pj Studying the terminal energy En or £'_„ (which 



is always an upper bound of the average final energy Ef) from Eqs. (81), (82 1, one sees that its maximum value is 
En ~ 20.1 GeV. This is however not enough to explain the arrival of much higher energetic neutrinos, so that one 
would still need to consider a matrix element with much stronger suppression than those considered in this work. 



Similarly, analyzing the decay length i in units of the distance L between CERN and Gran Sasso, we get from Eqs. ( 78 ) 



and (791 



i^K^_ A^ + nf ^17 GeV V" __^_ ^,,^^ _ f \p\ 



1 ,when e(b1)=(^l (94) 



'^ Assuming electron antineutrinos propagate with the same speed as the muon neutrinos detected by OPERA. 

* Hereafter in this section we will consider the second example of matrix element which is the most favorable for the superluminal 
interpretation of the OPERA results. 
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FIG. 2. Value of the terminal energy Eq (approximate final energy of a neutrino with initial energy higher than Eq) for different 
distances of propagation, as a function of the degree of superluminality in the neutrino speed. 



and 



where 



£_ 
L 



[l-af 



^- 



17GeV\ 



-3q 



E 



when e{\p\) ~ 



1.7 X 10-5 



[i;(17GeV)-lf V~^/ 



km 
IT 



(95) 



(96) 



and w(17GeV) — 1 = 2.35 x 10^^ is the velocity at 17 GeV required for the superluminal interpretation of the results 
reported by the OPERA collaboration. Eqs. (94) and (95) tell us that the decay length is much shorter than L when 
the energy is significantly larger than 20 GeV. This fact can be seen from the dimensionless factor r which is very 
small for these energies. For example, for E = hQ GeV, that is roughly a half of the maximum energy detected by 
OPERA, r = 5.7 x 10"'^ and the maximum value (with respect to n and a) of the decay length at this energy is 
/ = 2.8 X 10^^ • L, much shorter than the distance between CERN and Gran Sasso. 

The present analysis gives us an idea of the difficulties one finds to accommodate in a theoretical framework the 
results announced recently by the OPERA collaboration. In fact, the absence of energy loss in neutrino propagation 
is a more sensitive way to detect superluminal speeds. Taking our result for the second matrix element in the case 
of a momentum independent speed of propagation, the observation of the arrival of neutrinos of energies higher than 
60 GeV after a propagation of L = 730 km, gives a bound on the possible value of r\Q through Eq. ( 92 ) , 



riQ < 5.8 X 10" 



(97) 



which, according to Eq. (86), corresponds to a difference in the time arrival of neutrinos with respect to luminal 
speeds of Ato < 7.0 ns. This means that the non-observation of pair production in OPERA is incompatible (under 
the assumptions we are considering) with its time-of-flight measurement of ~ 60 ns, and that this non-observation is 
sensitive to tiny time differences, beyond the present precision of OPERA time measurements. 

This conclusion is of interest for any experiment trying to repeat the OPERA measurements. Superluminal speeds 
can be better detected through a deformation in the spectrum (of course, a direct measurement of the time-of-flight is 
also advisable, since it does not contain any assumption on the pair production process). In order to appreciate this 
deformation, we need to send neutrinos with energies higher than Eg, which depends on the distance of propagation 
and the degree of superluminality through Eq. ( 92 ) (we consider the case of a momentum independent speed as an 



illustrative case) . Therefore, an hypothetical experiment should consider both the neutrino energies of the beam and 
the distance of propagation in the sensitivity to the degree of superluminality. This sensitivity is shown in Fig. [2] in 
terms of {v — 1), the difference between neutrino and photon speeds, using our calculations with the second matrix 
element (the better motivated theoretically). Fig.^shows that, in order to be sensitive to {v— 1) of the order of 10"^, 
we should use neutrino beams of energies higher than 110 GeV if L ~ 730 kmr] while it would suffice to use neutrino 



In 2007 the MINOS collaboration reported a superluminal propagation speed v — 1 = (5.1 it 2.9) X 10 ^ (at 68% C.L.) for a (mostly 
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beams of energies higher than 70 GeV if L ~ 7000 km (the distance between CERN and the Soudan underground 
mine in USA). 

VII. CONCLUDING REMARKS 

In this work we have investigated the dependence of the charged lepton pair emission by superluminal neutrinos 
Vi — >■ Vi e~ e^ on the dispersion relation for neutrinos, and on the dynamical matrix element of the process. General 
expressions for an arbitrary dispersion relation and for various examples of matrix elements have been obtained. On 
the one hand, for a given dispersion relation, different choices of the matrix element lead to decay rates which differ 
by factors of order one. On the other hand, for a given matrix element, different choices of the dispersion relation 
lead to energy- independent, dimensionless factors (i.e., pure numbers) in the decay rates which can differ by one order 
of magnitude. These are the main new results of the present investigation. Estimates of the sensitivity of different 
observations of high energy neutrinos to a possible departure from SR kinematics in the neutrino sector have been 
presented. The dependence of the results on details of the theory (choice of modified dispersion relation and modified 
matrix clement) is an indication that high energy neutrino physics can be a very good laboratory to explore possible 
deviations from SR. 

Concerning the generality of the analysis, the assumptions on which it is based were already pointed out in the 
introduction: i) Rotational symmetry is preserved. Nevertheless, some works have built theories or models in which 
rotational invariance is not exact (see, for instance, Ref. [36]). ii) Energy and momentum are conserved in the conven- 
tional, additive way. However, during the last decade, there have been investigations which suggest the possibility of 
modifying these laws. Most of these explorations come from the quantum space-time realm and they deal with either 
deformations (see Ref. [SZ] for a review) or violations ([3S]) of space-time symmetries, iii) The relevant propagation 
speed of superluminal particles is the group velocity of their wave packets. Once more, several works have studied 
other alternatives both in canonical and in non-canonical space-times ( |39fHT| ) . 

While experimental results cannot be confirmed or refuted by theoretical investigations but by new experiments, 
what this present work shows is that the most likely possibility with respect to the observation reported by the OPERA 
collaboration is that either the observation of superluminal neutrinos is not confirmed or some of the assumptions 
indicated in the present work are not valid. In particular, in Doubly Special Relativity scenarios ([37'), where the 
energy-momentum conservation law is modified, forbidden processes in SR are generically forbidden too, as it was 
firstly pointed out in connection with the OPERA results in Refs. [I9 l l42 ] . 

Some works have studied the pair production decay width before us. In Refs. [T51I33], the decay width is computed 
for the case of constant velocity without specifying, however, what matrix element or relevant Lagrangian (if any) has 
been used to obtain the result. We reproduce that result for the first example of matrix element which does not come 
from a SU{2) gauge invariant underlying field theory. The pair production decay width is also computed in Ref. [22] 
for the momentum independent velocity case with two choices for a four-fermion interaction Lagrangian which are in 
fact equivalent to our second and fourth matrix elements. 

Despite the fact that we have concentrated on the pair production reaction, some of the techniques developed in the 
present work, and in particular the collinear approximation, can be extended to other processes like neutrino splitting 
^i ~^ ^z ^j i^j- In the case of a momentum dependent speed for neutrinos, splitting is a relevant mechanism for energy 
loss in the propagation of superluminal neutrinos. Other processes contributing to the energy loss of a propagating 
superluminal neutrino, such as the production of heavier charged leptons and the contribution of virtual W^ can also 
be treated within the same approximation. The results for these processes, as well as an analysis of the uncertainties 
in the calculation of the production of very high energy superluminal neutrinos with a general dispersion relation will 
be presented elsewhere [33]. 

NOTE ADDED 

While this manuscript was being completed, the OPERA collaboration announced the identification of two sources 
of error in the determination of the time of flight of neutrinos from CERN to Gran Sasso. One of them was a faulty 
connection in the optical fiber cable that brings the external GPS signal to the experiment master clock. That this 
was most probably the origin of the apparent neutrino superluminality was later confirmed (after this paper had been 
submitted for publication) by the results of muon measurements presented by LVD f45l , which show that the timing 



rtiuon) neutrino beam with an average energy of 3 GeV propagating 735 km between Fermilab and the Soudan mine j35| . It is expected 
that the MINOS team will report new results at higher precision in the near future. The OPERA, ICARUS, BOREXINO and LVD 
experiments at the Gran Sasso underground laboratory will also provide a more precise measurement of the neutrino speed. 
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between LVD and OPERA became misaligned from the middle of 2008, around the time that OPERA- 1 began, 
and remained stably misaligned by about 73 nanoseconds until the end of 2011, which is when the fiber problem 
was identified and eliminated. Almost at the same time, the ICARUS collaboration presented new results from the 
October-November 2011 campaign of measurements, which allowed a very accurate time-of- flight measurement of 
neutrinos from CERN to LNGS on an event-to-event basis, collecting seven neutrino events which are compatible 
with luminal speed |46j . Although new measurements are planned during 2012, the new information indicates that 
most probably the superluminal signal will go away, in agreement with the conclusions of the analysis presented in 
Section 6. 



ACKNOWLEDGMENTS 

This work is supported by CICYT (grant FPA2009-09638) and DGIID-DGA (grant 2010-E24/2). 



[1 
[2 
[3 

[4; 

[5 
[6 

[7: 

[9 
[10 

[11 

[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 

[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 
[33 
[34 
[35 
[36 
[37 
[38 

[39 
[40 
[41 
[42 



C. P. Burgess, J. M. Cline, E. Filotas, J. Matias, G. D. Moore, JHEP 0203, 043 (2002). [hep-ph/0201082]. 

V. A. Kostelecky, S. Samuel, Phys. Rev. D39, 683 (1989). 

J. Alfaro, H. A. Morales- Tecotl, L. F. Urrutia, Phys. Rev. Lett. 84, 2318-2321 (2000). [gr-qc/9909079] . 

F. Girelli, S. Liberati, R. Percacci, C. Rahmede, Class. Quant. Grav. 24, 3995-4008 (2007). [gr-qc/0607030] . 

C. Barcelo, S. Liberati, M. Visser, Living Rev. Rel. 8, 12 (2005). [gr-qc/0505065]. 

F. Girelli, E. R. Livine and D. Oriti, Nucl. Phys. B 708, 411 (2005) [gr-qc/0406100]. 

L. SmoUn, arXiv:0808.3765 [hep-th]. 

R. Aloisio, P. Blasi, P. L. Ghia and A. F. Grillo, Phys. Rev. D 62, 053010 (2000). [astro-ph/0001258]. 

D. Horns and M. Meyer, JCAP 1202, 033 (2012) [arXiv:1201.4711 [astro-ph.CO]]. 
U. Jacob and T. Piran, Phys. Rev. D 78, 124010 (2008) [arXiv:0810.1318 [astro-ph]]. 
T. Adam et al. [OPERA Collaboration], arXiv: 1109.4897 [hep-ex]. 

J. Ellis, N. Harries, A. MeregagUa, A. Rubbia and A. S. Sakharov, Phys. Rev. D 78, 033013 (2008). [arXiv:0805.0253 

[hep-ph]]. 

A. Chodos, A. L Hauser and V. A. Kostelecky, Phys. Lett. B 150, 431 (1985). 

E. Giannetto, G. D. Maccarrone, R. Mignani and E. Recami, Phys. Lett. B 178, 115 (1986). 

N. Arkani-Hamed, S. Dimopoulos, G Dvah and J. March-Russell, Phys. Rev. D 65, 024032 (2001) [hep-ph/9811448]. 
M. M. Anber and J. F. Donoghue, arXiv:1110.0132 [hep-ph]. 

D. M. Mattingly, L. Maccione, M. Galaverni, S. Liberati and G. Sigl, JCAP 1002, 007 (2010) [arXiv:0911.0521 [hep-ph]]. 
A. G. Cohen, S. L. Glashow, Phys. Rev. Lett. 107 (2011) 181803. [arXiv:1109.6562 [hep-ph]]. 

J. M. Carmona and J. L. Cortes, [hep-ph/ 11 10.0430]. 

L. Maccione, S. Liberati and D. M. Mattingly, arXiv:1110.0783 [hep-ph]. 

E. Ciufloh, J. Evslin, X. Bi and X. Zhang, arXiv:1112.3551 [hep-ph]. 

F. Bezrukov and H. M. Lee, Phys. Rev. D 85, 031901 (2012) [arXiv:1112.1299 [hep-ph]]. 

S. M. Carroll, J. A. Harvey V. A. Kostelecky C. D. Lane, T. Okamoto, Phys. Rev. Lett. 87, 141601 (2001). [hep- 

th/0105082]. 

J. Kowalski-Glikman, Phys. Lett. B 547, 291 (2002) [hep-th/0207279] . 

L. Susskind, J. Math. Phys. 36, 6377 (1995) [arXiv:hep-th/9409089] 

A. Shomer, arXiv:0709.3555 [hep-th]. 

G. F. Giudice, S. Sibiryakov and A. Strumia, arXiv:1109.5682 [hep-ph]. 

K. Hirata et al. [KAMIOKANDE-II CoUaboration] , Phys. Rev. Lett. 58, 1490 (1987). 

R. M. Bionta et al, Phys. Rev. Lett. 58, 1494 (1987). 

M. J. Longo, Phys. Rev. D 36, 3276 (1987). 

L. Stodolsky Phys. Lett. B 201, 353 (1988). 

R. Cowsik, S. Nussinov and U. Sarkar, Phys. Rev. Lett. 107, 251801 (2011) [arXiv:1110.0241 [hep-ph]]. 

X. -J. Bi, P. -F. Yin, Z. -H. Yu and Q. Yuan, Phys. Rev. Lett. 107, 241802 (2011) [arXiv:1109.6667 [hep-ph]]. 

L. Gonzalez-Mestres, arXiv:1109.6630 [physics. gen-ph]. 

P. Adamson et al. [MINOS CoUaboration], Phys. Rev. D 76, 072005 (2007) [arXiv:0706.0437 [hep-ex]]. 

A. G. Cohen and S. L. Glashow, Phys. Rev. Lett. 97, 021601 (2006) [arXiv:hep-ph/0601236]. 

G. Amelino-Camelia, Symmetry 2, 230 (2010) [arXiv: 1003.3942 [gr-qc]]. 

F. Mercati, D. Mazon, G. Amelino-Camelia, J. M. Carmona, J. L. Cortes, J. Indurain, C. Laemmerzahl and G. M. Tino, 

Class. Quant. Grav. 27, 215003 (2010) [arXiv: 1004.0847 [gr-qc]]. 

S. Mignemi, Phys. Lett. A 316, 173 (2003) [arXiv:hep-th/0302065]. 

M. Daszkiewicz, K. Imilkowska and J. Kowalski-Glikman, Phys. Lett. A 323, 345 (2004) [arXiv:hep-th/0304027]. 

S. Ghosh and P. Pal, Phys. Rev. D 75, 105021 (2007) [arXiv:hep-th/0702159]. 

G Amelino-Camelia, L. Freidel, J. Kowalski-Glikman and L. Smolin, arXiv: 11 10.0521 [hep-ph]. 



[43] Y. Huo, T. Li, Y. Liao, D. V. Nanopoulos and Y. Qi, Phys. Rev. D 85, 034022 (2012) [arXiv:1112.0264 [hep-ph]]. 

[44] J.M. Carmona, J.L. Cortes and D. Mazon, in preparation. 

[45] See A. Zichichi and G. Sirri talks in the mini-workshop "LNGS results on the neutrino velocity topic" 

http://agenda.infn. it/conferenceDisplay.py?confId=4896. 
[46] M. Antonello et al. [ICARUS Collaboration], arXiv: 1203.3433 [hep-ex]. 



19 



